Parameterised free functions and saddle stresses in modified gravity 
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Building on previous work, we explore the parameter space of free functions in non-relativistic 
modified gravity theories more widely, showing that in fact the two broad regimes present have 
similar functional forms between different models. Using different parameterisations, we investigate 
the effects on scaling tidal stresses as well as attempt to constrain the (hitherto poorly understood) 
deep MONDian scaling C. We also consider a new intermediate MOND limit in these theories and 
what it tells us about the transition between these regimes. Finally we suggest a model independent 
framework, with the aim of constraining the MONDian parameter space using future data, such as 
the forthcoming LISA Pathfinder mission. 



I. INTRODUCTION 



The concordance model of modern cosmology rests 
soundly on two cornerstones, a universe filled mostly with 
cold dark matter and dark energy (described by a cosmo- 
logical constant), ie ACDM, with underlying dynamics 
characterised by Einstein's theory of General Relativity 
(GR) . Whilst this model explains the early universe with 
every increasing accuracy, as long as there remains the 
lack of direct detection of a dark matter particle (baring 
unviable candidates such as neutrinos 1] ) , it remains pru- 
dent to consider alternatives. One such pathway avail- 
able is to modify the underlying dynamics themselves, 
subject to the condition that above certain scales we re- 
store our familiar Newtonian limit. MOdified Newtonian 
Dynamics [1, Q (MOND) provides just such a scheme. 
In the past decade, the potential accomplishments of 
MONDian theories have been put on an equal pedestal 
to GR with the development of fully relativistic mod- 
ified gravity theories, examples of such include Beken- 
stein's ground breaking theory of TeVeS [3], Einstein 
iEther theories Q and their generalised friends 
well as bimetric theories [l^ and various others [11 - 
Whilst the MONDian paradigm provides a useful frame- 
work to making connection to observables, the free func- 
tions and parameters in these theories remain relatively 
unconstrained, leading to a problem of fine tuning. Much 
work has been done investigating MONDian effects on 
the largest scales, for instance applying constraints from 
galactic data when seeking dark matter alternatives ^14- 
llSj . the much hailed Bullet cluster has been considered 
for what it can tell us about the necessity or needless- 
ness of dark matter and MOND [19- 23], as has applying 
Lorentz violating mechanisms (typically well constrained 
in the matter sector) to the gravity secto r [23 . [25[ . De- 
viations from the inverse square law have also 
been considered, however little more seems to be known 
about constraining modified gravity theories purely in 
the solar system. A chance of extending the forthcom- 
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ing LISA Pathfinder (LPF) mission [2i,[3g, to include 
probing the low acceleration regime around gravitational 
saddle points, provides just such an opportunity, both for 
testing and also cleanly constraining these theories. 

Previous work has explored the potential of using LPF 
measurements to test MONDian theories IsH]. We in- 
vestigated the expected signal to noise ratio (SNR) from 
LPF for such a test, as well as considering effects from 
a variety of instrumental noise models, trajectories past 
the saddle and issues arising from systematics such as 
self gravity [sl]. We also considered the possibility of 
constraining MOND in the event of a null result, as well 
as rescaling parameters in theories with preferred acceler- 
ation scale [33| . By investigating the symmetries present 
in such a scheme, we developed a simple algorithm for 
scaling tidal stresses in the event of parameters ao and 
K assuming values different those originally considered 
(short cutting a lot of computational work). In doing so, 
we were able to show broadly, that the scaling of parame- 
ters (within an order of magnitude) could potentially save 
or slaughter a dubious result. The effects from any viable 
cosmology [l^, 113] in these theories then becomes rele- 
vant, an example being applying BBN constraints to limit 
the effects of Gat renormalisation - which we can convert 
into a bound on variations to k. Similarly constraints 
from fifth force experiments and bounds on anomalous 
accelerations on solar system scales could feed into our 
theories, but these would appear as parameters in any 
viable free function - something beyond the scope of our 
simple scaling argument. 

In this work, we aim to investigate the effects of varying 
the free functions in these theories (in what have become 
known as Type I theories of MOND, see [12]). Here we 
show that a similar tidal stresses scaling algorithm exists 
in the linear regime of the theory, allowing us to rescale 
results previously computed from a particular (fiducial) 
model of /i and in fact further symmetries present here 
can be exploited for computational gain also. We also ex- 
plore the behaviour of the (somewhat mysterious) deep 
MONDian scaling C, which has been determined numeri- 
cally, but remains poorly constrained otherwise - here we 
consider how it gets rescaled between different models. 

The structure of this paper is as follows, firstly we con- 
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sider analytical solutions in the linear regime of the the- 
ory for a simple free function parameterisation, showing 
how the functional forms of solutions vary in different 
regimes and develop a new intermediate regime. In sec- 
tion mil we consider an application of this work to tran- 
sients for some more general fi and make considerations 
towards rescaling tidal stresses. We also attempt to con- 
strain the scalings in these models, comparing numerics 
in different models and providing some understanding 
as their behaviour. In the process, we develop a model 
independent framework for different free functions as a 
first step towards understanding the MONDian param- 
eter space. We finish with some broad conclusions and 
look towards future work. 



We also find the extra acceleration felt by test particles 
cannot be written down as a closed form expression in 
t/, however it can be given as 

.F = -V^.l^u(l + 1) (9) 

which in the quasi-Newtonian (QN) limit means « U" 
and in the deep MONDian (DM) hmit is z" « [/"/(n+i). 
With these results in mind, lets proceed to the same an- 
alytical solutions as before, using the linear Newtonian 
approximation. 

A. Quasi-Newtonian Regime 



II. ANALYTICAL RESULTS 

In Type I MONDian theories, dynamics in the non- 
relativistic regime result from the joint action of the usual 
Newtonian potential (associated with the metric) and 
a "fifth force" scalar field, 4>. The total potential acting 
on non-relativistic particles is their sum $ = $jv + 0- 
The MONDian field cj) is ruled by a non-linear Poisson 
equation: 



V • (^(z)V0) = kGp, 



(1) 



with z — -p-^-^^, here k is a dimensionless constant and 

47r ao ' 

ao is the usual MONDian acceleration. We begin from 
the linear MONDian field 



K V0 

U = M 

An ao 



(2) 



and will make use of the fact that since U — iiz, we can 
write any transition function also solely as /i = IJ-{U). We 
parameterise our yit function as 



(3) 



1 + z" 

and so we can write dimensionless vacuum equations 

V-U = (4) 
AmU^V AV + IJ A WU^ ^ (5) 



with 



ainC/2 2 

Am = — = - I n 

a In /i n 



1 



(6) 



where Am is picked for notational clarity later. Other pa- 
rameterisations will be considered later, but this choice 
represents a first step towards a more general free func- 
tion. Although we cannot write a closed form for m in 
terms of U, we can, as before, write down its form in the 
different regimes 

z « 1 ^l^z^^Am^ 2{n + l) 



z>l n^l ^ Am 



2C/" 



n 



(8) 



As before, we can separate out the QN field into a 
linear Newtonian and a curl like component which satisfy 



U 



Un + Uo 



Uo = -N(^) 



V • U2 



VAU, = - 



UqAVIUqP 
2|Uo|"+V" 



(10) 
(11) 
(12) 

(13) 



In spherical polar coordinates, we have the linear New- 
tonian profile 



N = A^,.e,. + iV^e^ 

Nr = i(l + 3cos2V') 
3 

~ — - sin 2-ip 
and using the notation 

U2 — Ur^r + U^e^ 

equations ([T^ and ([T^ take the form 

d 



19,0, 1 
- — {r^Ur) + — , „ , 
r or r sm o%p 



d_ 

dr 



dUr 







with 



-3n 



23"/2sin2V' 
(5 + 3cos2V')i+"/2 



(14) 
(15) 

(16) 
(17) 

(18) 
(19) 

(20) 



and so we have to discuss the effect of varying n on this 
extra curl force. 



1. The n = 1 case 



Here the source function becomes 
6\/2sin2'0 



Sl 



(5 -I- 3 cos 2^/^)3/2 



(21) 
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and equation (IT9|) suggests as ansatz of the form 

U2 = Bi(V') - (Fi(V')e, + Gi(V)ev,) (22) 

This reduces equation (IT^ to 

Gi = F[ + si (23) 

and hence equation (fT8)) to 

2Fi + cot V' + F'^ = - (s'l + si cot V') (24) 

We can solve this using the standard techniques of in- 
homogeneous ordinary differential equations (ODEs), to 
find expansions of F and G 

Fi « -0.23218- 0.7201 cos2?/' + 0.1306 cos4?/' 
Gi « 0.5115 sin2'(/'- 0.0556 sin 41/^ (25) 

We find the extra acceleration felt by test particles is 
hence given by 



SF 
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which for U ^ 1 suggests 



SF: 



U 



(26) 



(27) 



The first term, we call i5Fo, is simply a fully Newtonian 
term which renormalises the gravitation constant. The 
second term, which we denote <5Fi, is here just a rescaled 
unit vector of the Newtonian potential 



dFi 



167rao 



N(^) 



K V10 + 6COS 2-0 



(28) 



The third term, which we denote SF2, is just the curl 
field contribution 



2. The n = 2 case 



(29) 
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FIG. 1: The angular profile functions F and G giving the 
direction of the curl field in the QN region in a few 

different cases. We omit Gs here due to it vanishing for all -0. 



In this case, the source reduces to 

48 sin2?A 



Si 



(5 + 3 cos 20)2 



(30) 



which is just a rescaled solution of those from [351 • The 
form of equation (IT9t here suggests both C/^ and be- 
have as 1/r, hence we can rewrite our ansatz as 



B 



,(^) ^ (r^) (F2(0)e, + G2(0)e^) (31) 

Summarising the results of the calculation here, we find 
the ansatz collapses equation (jl9|) . allowing a simple sep- 
aration of the components of U: 



■A 



Gi sin -0 



5 + 3 cos 20^ 
tan-i(V3- 2tan|) -Htan-i(V3-H2tan|) 



with 



A = B 



in 



We can then insert this expression into the expansion in 
equation ([SJ for SF with [/ ^ 1 as before. 



3. The n = 3 case 



Here the source function becomes 
144\/2sin20 



S3 



(5 -I- 3 cos 20)5/2 
We find solutions using a separable ansatz of the form 



(33) 



V3 



+A cos il; + B 



(32) 



(^)'b3(0)= (^)'(F3(0)e, + G3(0)e^) (34) 
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with profile functions 

F3 = 



I6V2 



(5 + 3 cos 2^/^)3/2 



G3 = (35) 
which satisfy the boundary conditions. 

4. Other Cases 

For some power n > 1, we can make a more general 
ansatz as 



U 



2 — 



n-l 



(36) 



which we find now does not collapse equation ((T9)) . but 
rather by combining with (|18p. we get a second order 
sourced ODE for F„ in 

F„(n-2)(n-3) + i^>otV + K' = -« + s„cotV;) (37) 
and from (HH) 

(2 - n)Gn -Fl^^sn (38) 

The homogenous solutions of equation (j37p are simply 
Legendre polynomials in cos?/; of order {n — 2) and the 
full inhomogenous solution can be found using standard 
ODE techniques. We find a generic feature of solutions 
in this regime is that 



(5Fl+<5F2)(X-;l^ 



(39) 



and since further terms in (5F tail off ever faster, the 
relative importance of the curl terms becomes diminished 
in the large n limit. 



B. Deep MONDian Regime 

We can write equations @ and ^ here as 

d 



—7r{^ Ur) + — , „ , 

or r sm ip otp 



where a — 2 is used for notational convenience later and 
C is a constant required for matching between the two 
regimes. Applying this gives us a pair of coupled equa- 
tions for F and G 

G' + G coti^p) + aF = 0(45) 
pd{F^_jtSLl + 2[a'G-2mF']{F^ + G^) = 0(46) 



with 



a{2m - 1) + 2(1 - m) 



(47) 



where we should be clear to distinguish between the ra- 
dial exponent a and the extended variable a' . Since we 
are seeking solutions where the tidal stresses diverge, we 
can look at the additional force here 



SF = -S/cj) i 
which we rewrite as 



ro 



D 



(48) 



(49) 



where D is the angular profile in the DM regime. We see 
a < n + S yields divergent tidal stress solutions. Requir- 
ing n > 1 puts bounds on m as 



1 

— < TO < 1 



(50) 



by substituting in to and which after some manipulation, 
we find divergent solutions for 



a < 



6to- 2 
2to- 1 



and from inverting equation (|4^ 



a' + 2(to - 1) 
2to- 1 



(51) 



(52) 



4to / d{rUr 
r 



dU, 



r dtp ^ dr 



dr dtp 

which given the scaling symmetries of these equations 

U ^ U 

r — > A r 

suggests an ansatz for the potential as 

V^g{^] {F{^)Br + G[xP)e^) (44) 



this makes our bound a' < Am and hence given the 
bounds for to, we have that a' < 4 is always true. From 
similar considerations, we see that from requiring [/ 3> 1 
in the region r 3> rg, a > 2 in all cases (a point realised 
but not explicitly spelt out in [11]). These bounds are 
needed in picking out the particular a we require from the 
sequence which satisfy the equations and permit regular 
C/2=0(41) solutions. 



(sini/;C/^) = (40) 
d 



(42) 
(43) 



1. The m = 1 case 

If we consider solutions with to = 1 (equivalent to 
n = 1), then we are then guided to pick a' = a ~ 3.528 
(hence the results from ^35] stand) with profile functions 

Fi « 0.2442 + 0.7246cos2?/'-f 0.0472 cos 4?/' 

Gi « -0.8334sin2V'- 0.0368 sin 4^^ (53) 
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2. The i < m < 1 case 



Here, we have a' ^ a and so we need to find new 
solutions to equations ([45] -HH). We find solutions which 
neglect the derivative term in (j46|) as 



(54) 



F = acosip 
G = ^asmip 



(where a is a constant) here with values of a± given by 



1 - 2m 



(55) 
(56) 



but in fact the a_ solution only exists when m = 1, oth- 
erwise equations (|45] -l46 l) are not be simultaneously sat- 
isfied. We also find regular solutions exist for a discrete 
sequence of a(n)'s for each power n: {..., a_i, a±, ai, ...}, 
however now we have lifted the degeneracy that a_i = 
—ai (which only exists in the n = 1 case). 

For n = 2, we have solutions for a_i w —5.206, ai w 
3.983 and we pick solutions where 2 < a < 5, nicely 
selecting out ai, with angular profile functions as 

F2 w 0.248471 -I- 0.737261 cos2Vj-t- 0.05982 cos 4V' 
G2 ~ -0.9570sin2^- 0.057766 sin 41/^ (57) 

Similarly for n = 3, we have 2 < a < 6 and so ai ~ 
4.4057, with angular profile functions 

F3 w 0.257381 + 0.766674 cos2i/'-f 0.0624203 C0S4V' 
G3 « -1.1099sin2V'- 0.0756068sin4V' (58) 

and we compare a few angular profiles in Figure [H 



The m — ^ case 



In this large n limit, fi effectively becomes a step func- 
tion. We have a' — 1 and find equation (j46l) reduces 
to 



F 



d{F^ + 



2[G-F'](f2 + G2) = 



which can be further manipulated to just 

2 



1 



d_ fF_ 

di) \G 

with the simple solution of 

F = Gtan(?/'4-Gi) 



(59) 



(60) 



(61) 



and we pick Gi — ±7r/2 to satisfy the boundary condi- 
tions. Inserting this into equation (P5|) gives solutions of 
the form 



F 
G 



TG2C0sV'(sinV')" ^ 
G2(sini/;)"-i 




FIG. 2: A few numerically determined angular profile func- 
tions Fn and Gn in the DM region (dotted) compared with 
the Newtonian profile functions Nr and A*'^ (solid), in each 
figure respectively. Note the relative invariance of the radial 
profile and the slight changes in the azimuthal profile. 



This guides us to pick G2 = Tl and a = 2 here, giving 



F ~ cosTp 
G — ^psin-f/) 



(62) 



which mimics the a = 2 solutions seen previously. In this 
extreme case, the additional MONDian force is 



SF 



Anao D 



K D 

which gives tidal stresses of the form 



(63) 



(64) 



which obviously diverges as we approach the saddle (al- 
though the divergence is now relatively weaker than the 
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naive approach -"^ to MONDian tidal stresses suggests). 



C. An Intermediate Regime 



4- Type 11 Theories 

An important point to consider is what are the anal- 
ogous effects of changing free functions in other classes 
of MONDian theories. In the case of Type II theories, 
changes to each of the regimes are easy to include in the 
computation of the sourced Poisson equation [s^ . Given 
here 



47r 



V • [i^V$Ar] r Ci N • Vj/(r, Vj) (65) 



where v = v{v), v — |Uo| and Ci is some constant. We 
suggest solving this system with an ansatz 



(66) 



where C2 is another constant. We have a generalised DM 
limit of {F^)"-^^ (X Fjy, giving 

1 

V 

This fixes the radial exponent as 

n + 2 



n+1 



and also fixes the sourced ODE for the profile function 
a{a + 1)F„ + cot ^ F,', + F,'/ = g{^j, n) (67) 



where 



'l{ip,n) = (7 + 9 cos 2-0) 



(5 + 3cos2V')3"+2 



(68) 

Given the results of sourced ODEs like this from section 
III Al it seems unlikely therefore, that in this regime the 
solutions will remain similar, in stark contrast to the rel- 
ative invariance of the DM limit in the Type I theories. 
In the QN regime, a similar system of ODEs exist, with 
varying source functions and parameters depending on 
the form of the fall off from v ^ 1 and hence a simi- 
lar conclusion can be reached. A more detailed study of 
these generalised solutions however we leave for future 
work. 



^ The naive approach being to just take the 'rule of thumb' mod- 
ified inertia law, oc \/Tjv for a < ao, giving Sij oc r~'^ with 
e = i here, however as we have shown, for F^'^'^ oc Fjv with 
n > 1, we have solutions with 0.76 5; e < 1- 



Whilst we have a clear idea of the dynamics of U in 
the large and small acceleration regimes, we lack much 
detail in the 'near field' or intermediate regime (ie around 
z ~ 1 ) , except when we can estimate the size of the DM 
bubble. Our work thus far has focussed on finding the 
form and solutions to equation ([5]) in each limit, however 
really we only need to start at the form of the 4m and see 
how it scales in each limit. We can consider the leading 
order term to be of the form ~ If^ (where q is some power 
to be found) . We begin with the expression 



4to 



(69) 



where Cq is simply some dimensionless constant and first 
derive our results in the two well understood regimes. 
For the fi{n) models. 



2{n+l + z'' 
nCqZi 



1 



1 



and we seek solutions for ^ 
z regime. 



(70) 

1 in each limit. In the 



2(nH-l) 1 
nCq 



which can only approach unity when q = Q and 

_ 2(n+l) 
^'-^^ 

as before. Similarly in the z 3> 1 limit, 

2 



n-q 



hence for unity q = n and as before 



n 



(71) 



(72) 



(73) 



(74) 



The real power of this technique can be exploited to at- 
tempt to solve these models for z « 1, which has the 
expansion: 



(2 + ») ^,+i 
nCq 

2«(2n^ {2 + n)^q) 
nC„ 



(75) 



(z-1) 



+ 0{{z-iy) 
which we can solve at first order for g, 

2n 

q 



(2 + n)2 



and for unity we require 



Ca 



{2 + n) 



294 



(76) 



(77) 
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0.0 




FIG. 3: The ratio ^ plotted around z ^ 1 for different values 
of n. The blue (bold) line being n — 1, the purple (dashed) 
line n = 2, the red (solid) line n = 3 and the green (dotted) 
line n — 20. As we see, for small n, in this regime these 
models are relatively stable. 



Interestingly here, q is bounded for n > 1 as 



0<q< 



meaning that for large n, Cq — > 2, akin to the asymptotic 
DM regime here. We plot the full ^ profile to demon- 
strate the relative stability of this limit for different n's 
in Figure [31 Given however that here, 4m is not in gen- 
eral going to be constant, we are guided to pick QN- 
like perturbative solutions (at least for small n), taking 
U = Uo + U2 such that 



1-9 



B,(V') 



(78) 



consider ^'s of the form 



A* 



(1 -I- z'')'^/*' 
then as before, we compute 



(79) 



2 / 1 



a/b 



(80) 



where 



i U'' in the QN regime and z'' « [/VCf+i) in 
the DM regime. In the case of these free functions, our 
intermediate scaling parameter has the form 



2b 



q = 



{2 + af 



(81) 



which is clearly greater than unity for a = 1 and h > 
5, suggesting that the QN behaviour must have already 
been triggered before 2 ~ 1, since these intermediate 
solutions already display the asymptotic like behaviour 
(albeit with a 'stronger' curl term - in the z ^ 1 regime 
here, U2 (x r^^^). 

Our results are to be rescaled from our 'fiducial' model 
tJ'fiduciah presented in analytical and detailed numerical 
work previously [3l|, HI, . 



(82) 



or alternatively 



Vl + 4z4 - 1 

2^2 



(83) 



III. AN APPLICATION - TRANSIENTS 



One important use of this work is to look at how dif- 
ferent power law fall offs and transients from (p — -^^n 
effect our results. We can then examine how to rescale ex- 
isting templates for different /i functions. Whilst the /i(ri) 
model presents a nice parameterisation here, we should 
also consider multi-parameter families of functions, given 
that the various regimes are constrained by complemen- 
tary but different physical phenomena. Our 'usual' DM 
limit, fi on z, is motivated by the theory being a good dark 
matter replacement on low acceleration scales, a require- 
ment we can drop a priori, work with some generalised 
limit, II (X z"' , and then introduce later (as necessary). 
Whilst in the QN limit, the fall off from ^ — ^ 1 is gov- 
erned by agreement with fifth force experiments in the 
solar system (in particular the strong requirement of no 
anomalous accelerations on scales set by the orbit of Nep- 
tune [13 )• We see, therefore, it is prudent to consider at 
least a two parameter family of free functions. We can 



with 



4m 
SF 



4(1 + C/^) 

t/2 



1/4 



(84) 



In fact, if we pick some arbitrary model for fi{z), each 
regime takes the form 



z < 1 



1_V^ z»l 



(85) 
(86) 



n—q 



We then consider each a„ and 6„ as telling us the leading 
order (n — p, q) and higher terms of each expansion, with 
9 > 1 in all cases. In the DM regime, we find 



4m ~ 2 1 



p zP ap 



2(P+1) 



(87) 
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where p is the exponent of leading order term in fj.{z). 
We compute the term fj-fz^ from equation (|85l) . however 
at leading order this is obviously a constant - explaining 
why there is little variability in this limit. Similarly, in 
the QN limit, wc find (using U ~ z here) 




E 

n=q+l 



nbr, 



qbq 



where q is the exponent of leading order term in the ex- 
pansion of ^{z). As this result shows, here we are very 
much at the mercy of the free function we pick - naturally 
giving rise to the menagerie of solutions we found. 

Given these results, perhaps we can consider a simi- 
lar scaling argument for tidal stresses (as was considered 
in (33|). A naive approach is to consider a window func- 
tion (such as in section VC of [12]), ie one which preserves 
templates for r < tq (motivated by the relative invariance 
of the profile functions in this limit) but rescales them al- 
together differently outside of this. There are however a 
couple of issues that any rescaling algorithm needs to 
take into account: 

1. The Matching The constant C in the DM solu- 
tion is poorly constrained analytically and up to 
now has only been measured by looking at the ra- 
tio of numerical results to C = 1 analytical values 
(over the range r/ro = 0.05 — >■ 0.5). Care needs 
to be taken that the correct normalisation for any 
rescaled tidal stresses is found. 

2. Loss of Signal Since our analytical results are 
valid for r/rg <C 1, the nature of the window func- 
tion could cause us to loose signal when we rescale. 
Although these losses would be small compared to 
the signal deep inside the bubble, at the periphery 
and in the QN regime, noticeable losses could occur 
if a naive rescaling is done. 

To investigate such issues, we need to consider how the 
tidal stresses change when we vary the azimuthal and ra- 
dial components of the MONDian force (F0(r, -0)). Our 
anomalous tidal stresses are 

S,-^-J^ + ^^^ (89) 

however for simplicity here we will compute just Syy{x.), 
although with a suitable coordinate change, any tidal 
stress component could be picked. From the form of 
equation (|49p . we see we can write the MONDian force 
(in the linear regime) in the DM region as 



Fr 

F,i, 



where 



Ci r^g(V) 

a{n) - 2 

n + 1 
47rao 



(90) 



The tidal stresses therefore are 

Syy = ^r^-^(/(l + 7+(l-7)cos2^) (91) 



2 

+ (5(7-l) + /')sin27A-f 25'cos2^))-K^ 



where C2 = t^A is the rescaled Newtonian tidal stress at 

^ 47r 

the saddle. As Figure [2] shows, there is an approximate 
invariance of the profile functions in this regime, such 
that 

G„(0) ~ i{n)N^ (92) 

where ^(n) is some dimensionless linear scaling (for small 
n, we see ^ ~ 1). As this shows, any naive window func- 
tion rescaling of the tidal stresses runs the risk of ignoring 
relevant scalings dependent on "f{a{n)) and C. Similarly 
for the QN regime, in the linear regime (as seen in equa- 
tion (l?f)) ). the forces take the form 

Fr = C2iV,(V)r-HC3/(0)ri-" 
= C2iV^(0)r + C3 5(V')r-''" 

with 

_ 47rao 1 



(93) 



and hence the tidal stresses are 
C3 



Syy = ^r-"(/(2 + n(cosV^-l)) 

+ 2 cos ■(/'(sin ?/'(/' — ng) -I- cos?/;) 



(94) 



Here clearly we can try and play the same game with /, g 
- however we are hampered by the fact that the curl term 
has more variance between models. 

Another issue that must be addressed is the effect of 
the MONDian scaling C in these models. It would be 
telling if the scaling could indeed be a function of the 
parameters (perhaps n) in some way. Using the same 
techniques employed before, we reran our adaptive mesh 
code (presented in |3l|, but with simple adaptations, see 
Appendix]^ with smaller grids, to compute the ratio of 
analytical results [a prion with C = 1) and the numerical 
results for r/rp = 0.1 — >■ 0.5. As Figure 0] shows, the 
results are a little surprising. To try and make some sense 
of this, we need to look more closely at the matching 
between the regimes here. We find the existence of three 
types of relevant constant in our system of equations: 

• DM Regime For z <C 1, 4m — ^ C^*^, which be- 
comes relevant when computing the exponent a{n) 
in the DM regime solutions. 

• QN Regime For z > 1, 4m -J^ UP/C^^ where p 
is the leading order power relevant in the expansion 
of TO, such that when considering equation (jl3p . 



V AU2 



UqAVIUqI 
|Un|P+2 



VAUr 



(jQN 



(95) 
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FIG. 4: Deep MONDian scaling C for different ^ parameteri- 
sations, using different lattice sizes. We measured the ratio of 
the linear solutions compared to the full numerical solutions 
for r/ro = 0.1 — > 0.5 (since the smaller lattice sizes mean 
poorer resolution for small r, we raised the starting cutoff). 
As we see, whilst the errors reduce, the behaviour in n re- 
mains as we increase the lattice size (suggesting these effects 
are not artifacts). For comparison, we also present the origi- 
nal scaling for our fiducial /i, as well as a ^{a, b) model (with 
a = 1, 6 = 2). 



where U2 is the (renormalised) curl term and 
is the model dependent scaling. 

• Departures from Renormalised Gn 

Expanding /i in the z 3> 1 limit gives 



1 



z'P 



CI 



(96) 



where for consistency we match each by expand- 
ing out fj,~^ - we see that the coefRcients parame- 
terise the departures (at each order) from the usual 
Newtonian limit. 

We summarise the values of these for various parame- 
terisations: 



TABLE I: Scaling Constants for Various Models of fi 
C* fifiduciai fJ-{n) fJ-{a,b) 





4 


2{n + l) 


2(a + l) 




n 


a 




1 

4 


n 
2 


a 
2 




1 
4 


1 


a 
b 




1 

32 





a{a — b) 
26^ 




1 
128 





a(a — b){a — 2b) 
3b-i 



and we bring attention to the fact that in the ^(a, b) 
model, the DM regime parameter a plays a role in both 



regimes, whilst the fi{n) model is relatively constrained 
in the parameter space. Since we require the MONDian 
force be smooth and continuous, the matching between 
the different regimes must occur here. From equation ^ 
we see that 



SF = 



47rao U 

K /I 



which put together with these model independent pa- 
rameters, in the QN regime, reads as (the not unfamiliar 
expression of) 



U 



Cl(ri) 



Uo 



(97) 



where the higher order terms are 0{r^^'^P) or smaller and 
represent more complicated combinations of Uq and U2 
(see Appendix [Bl for more details). We matched this to 
the DM force 



— « Cf+i — 



D 

e 



(98) 



where H. is the leading order exponent in the expansion of 
jjL for z ^ 1. Here we choose the intermediate MONDian 
regime as 



U 



DM 



— ~ 0.5 ^ 1 

^0 



U, 



QN 



— ~ 1 4- 3 



Whilst matching the QN and DM regimes to find the 
scahng C might seem dubious, note that we seek not 
to predict the actual values of C{n) (since that would 
doubtless require knowledge of the system beyond linear 
order) but rather the scaling between different C's. In 
this way, we assume stricter validity in the linear regime 
but a more approximate one outside of this (as we have 
previously suggested with scaling rules). 

We summarise the results from such a matching and 
compare with numerical values: 



TABLE IL Comparison of Numerical and Matched C Scalings 



Cnum Cfnat ^num ^mat 



fJ' fiducial 

^(n = 1) 
^{n = 2) 
^(a = l;b = 2) 



0.8445 1.0523 

0.4042 0.5046 2.0893 2.0852 

0.8186 0.9847 1.0316 1.0686 

1.3163 1.7098 0.6416 0.6154 



where A# = {C fiducial /C^)^ represents the ratio be- 
tween the fiducial DM scaling and the value for each 
model of ^, in both the numerical and matched cases. As 
we see the values of C predicted in each case by match- 
ing typically overestimate its value. The ratio between 
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each value and the fiducial value (similarly obtaining by 
matching) is however much closer, suggesting we can find 
the scaling factor for the numerically obtained C fiducial 
by in general computing how it scales in this match- 
ing process. Naturally we should question that effec- 
tiveness of using these (so called) 'shooting' methods to 
interpolate between these disparate regimes. Such tech- 
niques are commonplace in the field of numerical relativ- 
ity, where matching asymptotic solutions to those close 
to a horizon, whilst maintaining regularity, is required. 
Recent work, however, has looked towards modelling cur- 
vature changes as akin to that of heat flows [H,[3i|. Such 
Ricci flow techniques allow for a much cleaner determi- 
nation of a systems dynamics, without the need for fine 
tuning of parameters. We leave the application of these 
spectral methods for future work. 

We can also consider free functions which do not 
asymptote to 1 (as explored in [l^, [13] )■ In this case 
we never enter the QN regime and so the DM rescalings 
described here will now apply to all the tidal stress. The 
effect on the scaling C however needs to be properly ex- 
plored, although naively we would expect it to take a 
value ~ 1, but we leave such issues to future work. 



IV. CONSTRAINING MOND 

A natural question to consider now is where do these 
results leave us when looking at data, such as measure- 
ments from LPF in the event of saddle flyby extension. 
How (if at all) can we constrain our theories from data? 
What exactly is the parameter space of MOND? 

The most likely scenario for an LPF test is a single sad- 
dle flyby, which given its (likely large) velocity, suggests 
data collected with be on the time scale of minutes. We 
posit that if any clear signal is seen above the Newtonian 
background, we can make the following inferences: 

• The 0*'* Order Approach, for 6 < rg, we assume 
that the signal is dominated by the DM regime, 
which given the relative invariance of the profile 
functions (assuming Fi{tp) ~ Fi, Gi{'tp) ~ Gi) 
means that the main scalings in the tidal stresses 
come from 7 and C. 

We can attempt to fit the DM tidal stresses to the 
signal by varying the value of 7 - remembering that 
really we have a{n) here and so really there only 
exists a strict series of 7(01(71)) for regular solutions 
of U. Once we find the correct 7, we can then 
consider the 'amplitude' of this signal, which can 
tell us (albeit broadly) about the QN regime, from 
the matching between (7^*^ — ?> Cf (ignoring C"^^ 
for the moment).^ 



For b > ro, the signal is sampling the QN regime, 
which means at lowest order 



U 



Un 



where 



h(^) = C', 



.m N(^) 
NP 



7-P-l 



(99) 



(100) 



Given that for small n, h{ip) ~ 0{1), we can first 
try fitting the radial fall off from the data and then 
once the exponent p is found, the various angular 
profile functions can be inferred. 



The 1** Order Approach, we introduce a cutoff 
in the signal, based on the impact parameter of 
the spacecraft (between the interior and exterior 
of the bubble). We fit the DM signal as before 
(now taking into account the scaling of Gi{ip)) and 
with the improved matching, make corrections to 
Cf and attempt to fit C^^ . 

We consider each parameter can be written as a 
perturbative expansion as 



+ . . 



(101) 



etc ... Allowing for improvements in parameter ac- 
curacy as we go up in approach. 



• The 2'"^ Order Approach, if we are blessed with 
plenty of tracking data and/or multiple flybys, we 
can make better determinations of our parameters. 

We note that not all of these parameters are independent, 
since a fixes both C^*^ and C^^, as well as the DM ex- 
ponent 7(a(n)). Similarly, constraining b fixes the forms 
of SFi cx Vq/Uq and SF2 oc U2(J7o,&). In this way, we 
can convert the distinctive MONDian signal of a positive 
result into a constraint on the MONDian parameter space 
akin to converting a negative result into a constraint on 
fi (as was considered in [32i]). 



^ Such a choice however is arbitrary, we can really choose either improve on its value and introduce the zeroth order term in the 

one of and to keep here and then at the next order discarded counterpart parameter. 
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V. CONCLUSIONS 

To conclude, we have presented a series of techniques 
for characterising, evaluating and rescaling the MON- 
Dian tidal stresses, that would be measured by LISA 
Pathfinder, should a saddle flyby be incorporated into 
the mission, for different models of /i function. Our goal 
was to detach these theories from their "alternative to 
dark matter duties" , considering two-regime functions 
with /X — >■ 1 at large 0, but /i oc z" (where in general 
n 1), when z is small. 

Broadly speaking, we find an approximate invariance 
for the angular profile functions in the DM limit irrespec- 
tive of the model of fi used, although changes do arise 
from the radial exponent 7 and DM scaling C. In the 
QN limit, we find a menagerie of solutions depending on 
the form of the fall off from /i — > 1. A very brief investi- 
gation into Type II MONDian theories suggests that the 
broad results here do not transfer over in to that class of 
theory and so any rescaling of tidal stresses must be done 
more carefully, in contrast to the invariance we showed 
in Type I theories. We also considered an intermediate 
MONDian limit, although it remains somewhat unclear 
what the solutions here can tell us about the transition 
between the other regimes. 

We suggested potential strategies to constrain the pa- 
rameter space of MONDian theories from data, identify- 
ing a framework of parameters that could be important 
for an experimental determination around saddle points. 
In doing so, we demonstrated the interplay between the 
DM and QN limits by using this framework to calcu- 
late how the DM scaling C varies in different models. 
We leave the computation of scalings in SNRs for future 
work, as well as other applications of these techniques. 



For completeness we compare with the expression for 

^fiducial • 



If K 

%2 " 2 



Appendix B: Extending the QN « DM Matching 
Calculation 



Since we need to match the MONDian forces in each 
regime (rather than say U), we start from the expression 

oc — 



along with 



and given that 



U = Vu-u 



l + 2CQ^^^^^-i^+(C«^) 



2 



C>(r-P) 




\ 1/2 
/ 

(Bl) 



we can put all of this together up to 3^ order as 
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Appendix A: Adaptations to the Numerical Code 

Minimal changes are required to run our code for these 
parameterised /i functions - the main difference being the 
term in the computation of the discrete divergence on 
the lattice (see Appendix A of [3l|), for the iJ,{n) models: 



dfj, n / K 



dg^ 2 V47rao 
and the /i(a, b) models: 



(Al) 



i9/i a 



u 



C>(ri) 



(B2) 



5(ri-P) 



Un -U?. 



r/2 



Uo 

^0 



(JQN 
•^0 



C>(ri-2p) 



■^0 



■^0 



U2 f |cfc«^ 



C>(ri-3P) 



^0 



TP 



C>(ri-3p) 



rP 



where p2 = p{p + 2)/2 and here the higher order terms 
are 0{r^^*P) or smaller and represent much more com- 
plicated combinations of Uq and 
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